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We discuss the structure of the "forward visible" (FW) parts of double and skewed distributions 
related to usual distributions through reduction relations. We use factorized models for double 
distributions (DDs) f(x,a) in which one factor coincides with the usual (forward) parton distri- 
bution and another specifies the profile characterizing the spread of the longitudinal momentum 
transfer. The model DDs are used to construct skewed parton distributions (SPDs). For small 
skewedness, the FW parts of SPDs H(x,£) can be obtained by averaging forward parton densi- 
ties fix — fa) with the weight p{a) coinciding with the profile function of the double distribution 
f{x, q) at small x. We show that if the x" moments /„ (q) of DDs have the asymptotic (1 — a^)"^^ 
profile, then the a-profile of f{x,a) for small x is completely determined by small- a; behavior of 
the usual parton distribution. We demonstrate that, for small f , the model with asymptotic pro- 
files for fn{a) is equivalent to that proposed recently by Shuvaev et al., in which the Gegenbauer 
moments of SPDs do not depend on f . We perform a numerical investigation of the evolution 
patterns of SPDs and gave interpretation of the results of these studies within the formalism of 
double distributions. 



I. INTRODUCTION 



Applications of perturbative QCD to deeply virtual Compton scattering (DVCS) and hard exclusive electropro- 
duction processes [Q-^ involve nonforward matrix elements (p — r\ 0(0, z) |p) | ^2^9 of quark and giuon light-cone 
operators. They can be parametrized by two basic types of nonperturbative functions. The double distributions 
(DDs) Fix,y;t) specify the Sudakov light -cone "plus" fractions xp'^ and yr+ of the initial hadron mo- 

mentum p and the momentum transfer r carried by the initial parton. Treating the proportionality coefficient 
C as an independent parameter one can introduce an alternative description in terms of the nonforward parton 
distributions (NFPDs) ^q{X; t) with X = x + y( being the total fraction of the initial hadron momentum taken by 
the initial parton. The shape of NFPDs explicitly depends on the parameter C characterizing the skewedness of the 
relevant nonforward matrix element. This parametrization of nonforward matrix elements by J-"^(X; t) is similar to 
that proposed by X. Ji who introduced off-forward parton distributions (OFPDs) (i, ^; t) in which the parton 
momenta and the skewedness parameter ^ = r^ /2P~^ are measured in units of the average hadron momentum 
P = {p + p')/2. OFPDs and NFPDs ||] can be treated as particular forms of skewed parton distributions 
(SPDs). One can also introduce the version of DDs ("a-DDs" Q) in which the active parton momentum is written 
in terms of symmetric variables k — xP -I- (1 -I- a)r /2. 

In our approach, DDs are primary objects producing SPDs after an appropriate integration. In refs. it was 
shown that using the the support and symmetry properties of DDs, one can easily establish important features of 
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SPDs such as nonanalyticity at border points X — (,0 [or x = ±^], polynomiahty of their and moments 
in skewedness parameters C and ^, etc. 

The physical interpretation of DDs F[x,y;t = 0) [or f(x,a;t =^ 0) ] and their relation to the usual parton 
densities f{x) suggests that the x-profile of DDs F{x,y), f{x,a) is driven by the shape of f{x) while their y or 
a-profile is analogous to the shape of two-body distribution amplitudes like [y(l — x — y)]", [(1 — x)"^ — a^]". Fixing 
the profile parameter n gives simple models [|^,^ for DDs which can be converted into models for SPDs. 

In the present paper, our main goal is to study the self-consistency of these models with respect to the pQCD 
evolution. In Section II, we briefly review the basic elements of the formalism of double distributions, discuss their 
support and symmetry properties and relation to usual parton densities. In Section III, we describe factorized 
profile models for double distributions and give explicit model expressions for skewed distributions. In Section IV 
we consider a practically important case when skewedness parameters C or ^ are small. The factorized models for 
DDs in this case can be taken in a very simple form f(x,a;t = 0) = f{x)p(a), where p{a) is a normalized profile 
function. As a result, SPDs T(^{X) (or H{x,£^)) in this model are obtained by averaging the relevant forward 
distribution f{x) over the interval {X — CiX) [or {x — £,,x + ^)] with the weight p{a) (we use the convention ^ 
that "tilded" parton distributions are those defined on the (—1,1) interval). In Section V, we study the impact 
of the pQCD evolution on the profile function p{a). Since the a-DDs are hybrids which look like usual parton 
densities wrt x and like distribution amplitudes wrt a, the simplest renormalization properties at one loop have 
the combined x"^Cf^'^~^"^{a) moments of f{x,a). As a result, independently of the initial condition, the a-profile 
of the x" moment /«(«) of the a-DD f{x,a) under the pQCD evolution asymptotically tends to (1 — a^)"+^. 
We investigate the "asymptotic profile model" in which /«(«) are given by their asymptotic form and show that 
it imposes a remarkable correlation f{x,a) = F{x/{1 — a^)) between the a;-dependence of the a-DDs and their 
a-profile. To study the impact of pQCD evolution on the DD based models of SPDs, we perform an explicit 
numerical evolution of SPDs. In Section VI, we describe a simple algorithm for the leading-log evolution of SPDs 
based on direct iterative convolutions of evolution kernels Wi^{X,Z) with SPDs !F(^{Z). In section VII, we discuss 
the results of our numerical calculations. In Appendix A, we show that the approximation (used in Ref. Q) in 
which the Gegenbauer moments of SPDs do not depend on skewedness, is equivalent to the asymptotic profile 
model for DDs. In Appendix B, we present explicit form of evolution equations for SPDs used in our numerical 
calculations. 



II. DOUBLE DISTRIBUTIONS 



In the pQCD factorization treatment of hard electroproduction processes, the nonperturbative information is 
accumulated in the nonforward matrix elements {p — r\ 0(0, z) \p) of light cone operators 0(0, z). For z^ — the 
matrix elements depend on the relative coordinate z through two Lorentz invariant variables (pz) and (rz). In 
the forward case, when r = 0, one obtains the usual quark helicity-averaged densities by Fourier transforming the 
relevant matrix element with respect to (pz) 

(p, s' I MO)mO, A)Mz) \p,s)\ ,2=0 = u{p. s')zu{P, s) (e-'^(^'^)/a(x) - e^"^^'^)/^^)) dx , (2.1) 

where E{0,z]A) is the gauge link, u{p' , s'),u{p, s) are the Dirac spinors and we use the notation jaz" = z. 
The functions fa{x), fa{x) can be also treated as components of the "tilded" distribution fa{x) = fa{x)9{x > 
0) — fa{x)9{x < 0), whose support extends to [— 1 < a; < 1]. In the nonforward case, we can use the double Fourier 
representation with respect to both (pz) and (rz) 

{p\ s' I MO)zE{0, z; A)M^) I | ,2=0 (2.2) 

= u{p', s')zu{p, s) dy j ^-i<P^)-^v'^^-) paix, y]t)e{0 <x + y < 1) dx + "0(r)-terms" , 

where the "0(r)-terms" stands for contributions which have the structure u{p' , s')r^^ z'^a^,^u{p, s)^{{pz), (rz)) [|l[ 
and {rz)u{p' , s')u{p, s)^'((rz)) [p3|p| and vanish in the r ^ limit. 
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For any Feynman diagram, the spectral constraints ~l < x < l,0<y< l,0<a; + y< 1 were proved in 
the a-representation |^ using the approach of Ref. jlOj. The support area for the double distribution Fa{x,y;t) is 
shown on Fig.|l|a. 
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FIG. 1. a) Support region and symmetry line y — x/2 for j/-DDs F{x,y;t). b) Support region for a-DDs f{x,a). 



Taking the r — hmit of Eq. (2.2), one obtains "reduction formulas" relating the double distribution Fa{x, y; t 
0) to the quark and antiquark parton densities 



Fa{x,y;t = 0)\x>o dy = faix) 



Faix,y;t = 0)\x<o dy = -/a(-a;) . 



(2.3) 



Hence, the positive-x and negative- a; components of the double distribution Fa{x, y; t) can be treated as nonforward 
generalizations of quark and antiquark densities, respectively. If we define the "untilded" DDs 



Faix, y; t) = Faix, y;t)\x>o ; Faix,y;t) = -Fai-x,l - y;t)\x<o , 

then X is always positive and the reduction formulas have the same form 

-X 

Fa,aix,y;t ^ Q)\x^(i dy = ,fa,aix) 



(2.4) 



(2.5) 



in both cases. The new distributions both "live" on the triangle Q < x,y < 1, < x + y < 1. Taking z in 
the lightcone "minus" direction, we arrive at the parton interpretation of functions Fa^aix,y',t) as probability 
amplitudes for an outgoing parton to carry the fractions xp^ and yr^ of the external momenta r and p. The 
double distributions Fix,y;t) are universal functions describing the flux of and r+ independently of the ratio 
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FIG. 2. a) Parton picture in terms of j/-DDs; b) meson-like contribution; c) Polyakov- Weiss contribution; d) parton 
picture in terms of a-DDs. 

The functions Fix, y; t) may have singular terms at x = proportional to Six) or its derivative(s). In such terms 
the sensitivity of the parton distribution to the plus-component of the initial hadron momentum is lost, and they 
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have no projection onto the usual parton densities. We will denote them by Fm{x, y; t) — they may be interpreted 
as coming from the t-channel meson-exchange type contributions (see Fig|2|3). As shown by Polyakov and Weiss 
, the terms with quartic pion vertex (Fig.ljc), in which the dependence on (pz) is also lost, correspond to an 
independent {rz)u{p' , s')u{p, s)\E'((rz)) type contribution. Both the meson-exchange and Polyakov- Weiss terms are 
invisible in the forward limit, hence the existing knowledge of the usual parton densities cannot be used to constrain 
these terms. In what follows, we consider only the "forward visible parts" of SPDs which are obtained by scanning 
the X parts of the relevant DDs. 

To make the description more symmetric with respect to the initial and final hadron momenta, we can treat 
nonforward matrix elements as functions of (Pz) and (rz), where P = {p + p')/ 2 is the average hadron momentum. 
The relevant double distributions fa{x, a ; t) [which we will call a-DDs to distinguish them from y-DDs F{x, y; t)] 
are defined by (l|J| 

{p'\M-z/2)zMzm ^u{p')zu{p) f dx C e-"(^^)-'"(''^)/2j=a(x,a;t)da +"0(r)~terms". (2.6) 

The support area for fa{x,a;t) is shown in Fig.^. Again, the usual forward densities fa{x) and fa{x) are given 
by integrating fa{x, a ; t = 0) over vertical lines x ~ const for x > and x < 0, respectively. Due to hermiticity 
and time-reversal invariance properties of nonforward matrix elements, the a-DDs are even functions of a: 

fa{x,a;t) = fa{x,-a;t). 

For our original y-DDs Fa,a{x,y;t), this corresponds to the "Munich" symmetry with respect to the interchange 
y ^ 1 — x — y established in Ref. . The a-quark contribution into the flavor-singlet operator can be parametrized 
either by y-DDs F^{x, y; t) or by a-DDs (x, a ; t). The latter are even functions of a and odd functions of x: 

f!{x, a; t) = {fa{\xl \a\-t) + h{\x\, \a\-t)} sign(x) -I- fi,{x, a; t) . (2.7) 

The valence quark functions fY {x, a ; t) are even functions of both a and x: 

f^{x,a;t) = U\xl\a\;t) - f-ai\xl\a\;t) + fl,{x,a;t) . (2.8) 

It is convenient to define the gluonic a-DD f'^ix, a; t) in such a way that its integral over a for t = 0, also gives 
the usual forward gluon density f^{x): 

{P - r/2 I z^z,G;:j-z/2)£;,b(-z/2, z/2; A)Gi,{z/2) \ P + r/2) \ ,2=0 (2.9) 

.1 ,1-1x1 

= u{p')zu{p){Pz) dx e-"^^''>-'°'^"'>/^ xf'^{x,a;t) da + 0{r) -terms. 

Jo 

The gluon SPD iJ'~^(5;, ^; /i) is constructed in this case from xf^{x,a;t). Just like the singlet quark distribution, 
the function f'~^{x, a;t) is an odd function of x. 

III. MODELS FOR DOUBLE AND SKEWED DISTRIBUTIONS 

The reduction formulas and interpretation of the x-variable as the fraction of p (or P) momentum suggest that 
the profile of F{x, y) (or f{x, a)) in x-direction is basically determined by the shape of f{x). On the other hand, the 
profile in y (or a) direction characterizes the spread of the parton momentum induced by the momentum transfer 
r. In particular, since the a-DDs /(x, a) are even functions of a, it make sense to write 

f{x,a) = h{x,a)f{x), (3.1) 

where h(x, a) is an even function of a normalized by 
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h(x,a)da — 1. (3-2) 

l+\x\ 

We may expect that the a-profile of h(x,a) is similar to that of a symmetric distribution amphtude (DA) ip{a). 
Since |a| < 1 — |a;|, to get a more complete analogy with DA's, it makes sense to rescale a as a = (1 — \x\)P 
introducing the variable (3 with x-independent limits: — 1 < /? < 1. The simplest model is to assume that the 
profile in the /^-direction is a universal function g{/3) for all x. Possible simple choices for g{f3) may be 5{P) 
(no spread in /3-direction) , |(1 — /3^) (characteristic shape for asymptotic limit of nonsinglet quark distribution 
amplitudes), y|(1 — /3^)^ (asymptotic shape of gluon distribution amplitudes), etc. In the variables x, a, this gives 

These models can be treated as specific cases of the general profile function 

/.(^)(x a) = r(26 + 2) [(l~N)^-aT .3 
^2;, a; 22''+ir2(6 + l) (1 - |a;|)2''+i ' 

whose width is governed by the parameter b. 
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FIG. 3. Integration lines for integrals relating SPDs and DDs. 

The coefficient of proportionality C — jp^ (or ^ = r+/2P+) between the plus components of the momentum 
transfer and initial (or average) momentum specifies the skewedness of the matrix elements. The characteristic 



feature implied by representations for double distributions [see, e.g., Eq.(2^)] is the absence of the ^-dependence 
in the DDs F{x, y) and ^-dependence in /(x, a). An alternative way to parametrize nonforward matrix elements of 
light-cone operators is to use C (or ^) and the total momentum fractions X = x-\-yC (or x = x + ^a) as independent 
variables. Integrating each particular double distribution over y gives the nonforward parton distributions 



j^i{x) ^ f dx f 

Jo Jo 



S{x + Cy~X)F,ix,y)dy 

rX/C rX/C 

= e{X > C) / F,{X- yC, y) dy + 9{X < C) / F,{X - yC, v) dy , (3.5) 
Jo Jo 

where C, = I — (■ The two components of NFPDs correspond to positive {X > () and negative {X < Q values of 
the fraction X' = X — ( associated with the "returning" parton. As explained in refs. [§,D, the second component 
can be interpreted as the probability amplitude for the initial hadron with momentum p to split into the final 



5 



hadron with momentum (1 — ()p and the two-parton state with total momentum r = C,p shared by the partons in 
fractions Yr and (1 — Y)r, where Y = X/C,. 

The relation between "untilded" NFPDs and DDs can be illustrated on the "DD-life triangle" defined by < 
x,y,x + y < 1 (see Fig.^). Specifically, to get T(;{X), one should integrate F{x,y) over y along a straight hue 
X = X — (^y. Fixing some value of C, one deals with a set of parallel lines intersecting the a;-axis at x = X . The 
upper limit of the y-integration is determined by intersection of this line either with the line x + y = 1 (this happens 
if X > C) or with the y-axis (if X < (). The line corresponding to X ^ ( separates the triangle into two parts 
generating the two components of the nonforward parton distribution. 

In a similar way, we can write the relation between OFPDs H{x,S,;t) and the a-DDs f{x,a;t) 

/•I /-i-l^l 

H{x,^;t) = / dx / S{x + — x) f{x,a;t) da . (3-6) 

J-l J-l + |a:| 

It should noted that OFPDs as defined by X. Ji |l[ correspond to parametrization of the nonforward matrix element 
by a Fourier integral with a single common exponential with — 1 < ai < 1, i.e., OFPDs are equivalent to "tilded" 



NFPDs. The delta-function in Eq.(3.6) specifies the line of integration in the {a;, a} plane. For definiteness, we 
will assume below that ^ is positive. 

Information contained in SPDs originates from two physically different sources: meson-exchange type contribu- 
tions J-q'^{X) coming from the singular a; = parts of DDs and the functions !F^{X), Tq{X) obtained by scanning 
the x 7^ parts of DDs F°-{x, y), F°'{x, y). The support of exchange contributions is restricted to < X < (. Up 
to rescaling, the function !Fq'\X) has the same shape for all C- For any nonvanishing X, these exchange terms 
become invisible in the forward limit C ^ 0. On the other hand, the support of functions JF?(X), !Fq{X) in general 
covers the whole < X < 1 region. Furthermore, the forward limit of such SPDs as Tq'°'{X) is rather well known 
from inclusive measurements. Hence, information contained in the usual (forward) densities /"(cc), /"(x) can be 
used to restrict the models for J^^{X), J^^{X). 

Let us consider SPDs constructed using simple models of DDs specified in Section III. In particular, the model 
/(°°)(a;,Q;) = 5{a)f{x) (equivalent to F'-°°^{x,y) = 6{y - x/2)f{x)), gives the simplest model H''°°\x,£,;t = 0) = 
f{x) in which OFPDs at t = have no ^-dependence. For NFPDs this gives 



i ^ ^ l-C/2 V l-C/2 

i.e., NFPDs for non-zero ( are obtained from the forward distribution f{X) = J-',^=o{X) by shift and rescaling. 

In case of the 6 = 1 and b = 2 models, simple analytic results can be obtained only for some explicit forms of 
f{x). For the "valence quark" -oriented ansatz /'^^)(a;,a), the following choice of a normalized distribution 

f^'\x)^^^-^x-^{l-xr (3.8) 
61 (1 — a) 

is both close to phenomenological quark distributions and produces a simple expression for the double distribution 
since the denominator (1 — x)^ factor in Eq. (O) is canceled. As a result, the integral in Eq. (3^) is easily 



performed and we get 

^y(5,e)l|5|>c = ^ (l - ^) { [(2 - a)ai - S:){xl'^ + xl'^) + [e 5)(a;?-" - xl'^)] e{i) + {i ^ -i)] (3.9) 



for |a;| > ^ and 



i/Mi,0l|5|<C = ^ (l - i) {^i""[(2 - + le x)] + (i ^ -i)} (3.10) 

in the middle — ^ < x < region. We use here the notation xi ~ (a; + ^)/(l + ^) and X2 = (a; — ^)/(l — ^) To 
extend these expressions onto negative values of ^, one should substitute ^ by |^|. One can check, however, that 
no odd powers of |^| would appear in the x^ moments of H^^ {x, ^). Furthermore, these expressions are explicitly 
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non-analytic for x = This is true even if a is integer. Discontinuity at a; = however, appears only in 
the second derivative of i.e., the model curves for H^^{x,^) look very smooth (see Fig.||). The explicit 

expressions for NFPDs in this model were given in ref. ||^. The relevant curves are also shown in Fig.^. 





0.2 0.4 0.6 0.8 1 -1 -0.5 0.5 1 

FIG. 4. Valence quark distributions; untilded NFPDs F^{x) (left) and OFPDs Hy{x,£,) (right) with a = 0.5 for several 
values of namely, 0.1, 0.2, 0.4, 0.6, 0.8 and corresponding values of ^ = C/(2 ^ C)- Lower curves correspond to larger 
values of 



For a = 0, the x > part of OFFD has the same x-dependence as its forward limit, differing from it by an 
overall (^-dependent factor only: 



(1-1^ 



|i|>0 +2 



(1 + 0= 



\i\<0- 



(3.11) 



The (1 



behaviour can be trivially continued into the |i| < ^ region. However, the actual behaviour of 



H^^ {x, £,)\a=o in this region is given by a different function. In other words, H^^ (x, ^)\a=o can be represented as a 
sum of a function analytic at border points and a contribution whose support is restricted by \x\ < ^. It should be 
emphasized that despite its DA-like appearance, this contribution should not be treated as an exchange-type term. 
It is generated by regular x part of DD, and, unlike <y9(a;/^)/^ functions changes its shape with ^ becoming 
very small for small 
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FIG. 5. Singlet quark distributions: tilded NFPDs (x) (left) and OFPDs H^{S:,0 (right) for several ( values 0.2, 
0.4, 0.6 and corresponding values oi ^ = (/{2 ~ (). Lower curves correspond to larger values of Forward distribution is 
modeled by (1 — x)"^ /x. 



For the singlet quark distribution, the a-DDs f^{x,a) should be odd functions of x. Still, we can use the 
model like ( p.8D for the x > part, but take /'^(x, a)|a;^o = ^ f''^H\x\: a) sign{x). Note, that the integral (3.6) 
producing H^{x,£,) in the |i| < C region would diverge for a ^ i/^ if a > 1, which is the usual case for standard 
parametrizations of singlet quark distributions for sufficiently large Q^. However, due to the antisymmetry of 
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f^{x,a) wrt X — > —X and its symmetry wrt a — > —a, the singularity at a = x/^ can be integrated using the 
principal value prescription which in this case produces the x - 
For a = 0, its middle part reduces to 



■X antisymmetric version of Eqs.(3.9) and (3.1C) 
3^2 _ 2x^^ ~ x^ 



(3.12) 



The shape of singlet SPDs in this model is shown in Fig. || 



IV. SPDS AT SMALL SKEWEDNESS 



To study the deviation of skewed distributions from their forward counterparts for small ^ (or Q) , let us consider 
the x> part of H{x, ^) [see Eq.(3.6)] and use its expansion in powers of ^ 



1 /-(i-^) d^f{i,a) 

(l-x) 



dx^ 



a' da +{1- if P^^"'"^-2-^^("'") 



V da 



dx 



a — l — x 



(4.1) 



where /(i) is the forward distribution. For small the corrections are formally 0(^^). However, if f{x,a) has a 
singular behavior like x^"", then 



d^f{S:,a) a{l + a) 
dx"^ x"^ 



fix, a) , 



and the relative suppression of the first correction is 0(^^/i^). Though the corrections are tiny for i 3> in the 
region a; ^ it has no parametric smallness. It is easy to write expicitly all the terms which are not suppressed in 
the 5: ~ ^ — ^ limit 



k=0 



{2ky. 7_i 952fe 



X — ^a, a) da + . . . , 



(4.2) 



where the ellipses denote the terms vanishing in this limit. This result can be directly obtained from Eq. ( |3.6| ) by 
noting that for small x, we can neglect the x-dependence in the limits ±(1 — |a;|) of the a-integration. Furthermore, 
for small x one can also neglect the x-dependence of the profile function h{x, a) in Eq. (^^) and take the model 
f{x,a) — f{x)p{a) with p{a) being a symmetric normalized weight function on —1 < a < 1- Hence, in the region 
where both x and ^ are small, we can approximate Eq. (3.6) by 



Hii; C) = "P" / fii ~ ^a)p{a) da + 



(4.3) 



i.e., the OFPD H{x; £,) is obtained in this case by averaging the usual (forward) parton density f{x) over the region 
X — £,<x<x + S, with the weight p{a). The principal value prescription "P" is only necessary in the case of 
singular quark singlet distributions which are odd in x. In terms of NFPDs, the relation is 



TciX) = "P" / fix - C(l + a)/2)p(a) da + 



(4.4) 



i.e., the average is taken over the region X ( < x < X . 

In fact, for small values of the skewedness parameters ^, C, Eqs. ( |4.3| ), (4.4) can be used for all values of x and 
X: if i; ^, Eq. ( ^ ) gives the correct result H{x;£_) = /(i) + 0(^2). Hence, to get the "forward visible" parts 
of SPDs at small skewedness, one only needs to know the shape of the normalized profile function p{a). 

The imaginary part of hard exclusive meson electroproduction amplitude is determined by the skewed distribu- 
tions at the border point. For this reason, the magnitude of [or H{£^, ^)] and its relation to the forward densities 
f{x) has a practical interest. This example also gives a possibility to study the sensititivity of the results to the 
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choice of the profile function. Assuming the infinitely narrow weight p{a) = 6{a), we have J-(^{X) — f{X — (/2) + . . . 
and H{x,S,) = f{x). Hence, both J^^iC) and H{^,£_) are given by f{xBj/2) since ( = XBj and ^ = + . . .. 

Since the argument of f{x) is twice smaller than in deep inelastic scattering, this results in an enhancement factor. 
In particular, if f{x) ~ for small x, the ratio Tl{C) = ^c(0/ fiO is 2". The use of a wider profile hmction 
p{a) produces further enhancement. For example, taking the normalized profile function 



and f{x)^x ° we get 



, , r(6 + 3/2) , r(26 + 2) . 2^fc .^ 

" r(i/2)r(5 + i) (^ - " ) = 25S^WTT) " " ^ (^-^^ 
^(.) ... = :1!(C) _ r(26 + 2)r(b-a + i) 

/(C) " r(26-a + 2)r(6+l) ^^-"^ 
which is larger than 2° for any finite b and < a < 2. The 2° enhancement appears as the n ^ oo limit of Eq.(E^). 



For small integer n, Eq.(45) reduces to simple formulas obtained in refs. |3,|7|. For n — 1^ we have 



— = = , (4.7) 

/(C) (l-a/2)(l-a/3)' ^ ^ 



which gives the factor of 3 for the enhancement if a = 1. For 6 = 2, the ratio (4.5) becomes 

^r^^(c) I 

/(C) (l-a/3)(l-a/4)(l-a/5)' 



(4.8) 



producing a smaller enhancement factor 5/2 for a = 1. Calculating the enhancement factors, one should remember 
that the gluon SPD !Fq{X) reduces to Xfg{X) in the C — limit. Hence, to get the enhancement factor corre- 
sponding to the fg{x) ^ x'^ small- a; behavior of the forward gluon density, one should take a = A — 1 in Eq.(p~^), 
i.e., despite the fact that the l/x behavior of the singlet quark distribution gives the factor of 3 for the TZ^^C) 
ratio, the same shape of the gluon distribution results in no enhancement. 

Due to evolution, the effective parameter a characterizing the small-x behavior of the forward distribution is an 
increasing function of Q^. As a result, for fixed h, the TZ^^\C,) ratio increases with Q^. In general, the profile of 
/(x, a) in the a-direction is also affected by the pQCD evolution. In particular, in ref. it was shown that if one 
takes an ansatz corresponding to an extremely asymmetric profile function p{o) 6(1 + a), the shift of the profile 
function to a more symmetric shape is clearly visible in the evolution of the relevant SPD. In the next section, we 
will study the interplay between evolution of x and a profiles of DDs. 



V. QCD EVOLUTION AND PROFILE OF DDS 

Both the shape of the forward distributions f{x; p) reflected in the x-dependence of the DDs /(cc, a; p) and 
their profile in the a-direction are affected by the pQCD evolution. At the one-loop level, the solution for QCD 



evolution equations is known in the operator form |18|, so that choosing specific matrix elements one can convert 
the universal solution into four (at least) evolution patterns: for usual parton densities {{p\ - ■ ■ \p) case), distribution 
amplitudes (0| . . . \p) case), skewed and double parton distributions {{p — r\ . . .\p) case). Since all the types of the 
pQCD evolution originate from the same source, one may expect an interplay between the x- and a- aspects of the 
DDs evolution. 

In the simplest case of flavor-nonsinglet (valence) functions, the multiplicatively renormalizable operators were 
originally foimd in Ref. ||l^ 

©.^■^ = (za+)" V^A"zC3/2(z B lzd+)i> . (5.1) 

In contrast, the usual operators ?/;A"z(z D)"'ip mix under renormalization with the lower spin operators 
{zd+)"-''^\''z{z D)''ip. The symbolic notation (z B / zd+) with B=B - B , 9+ =B + B=a + 5 and Cl''^{a) 
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being the Gegenbauer polynomials is borrowed from Ref. M]. In Ref. [LOl it was also noted that these operators 



coincide with the free-field conformal tensors. As pointed out in Ref. 
O^^ operators means that the Gegenbauer moments 



|, the multiplicative renormalizability of 

A^) = C j ^ Cl/\z/0 H'^'^iz, t A*) dz (5.2) 
of the skewed parton distribution H^^ {z,£,; fi) have a simple evolution j^: 



luA^o/Ar-/''" 



In fi/A 



(5.3) 



where /3o = 11 — ^Nf is the lowest coefficient of the QCD /3-function and 7„'s are the nonsinglet anomalous 
dimensions |2^,^. Going from SPDs to DDs, writing the SPD variable x in terms of DD variables x — x + a£, and 
using 

(-/^ + a)=± + (2./0"-' Cf/^^"-'(a) , (5.4) 

^ r(3/2)(n~ /)! 

one can express the Gegenbauer moments C„(^, /i) in terms of the combined [a;-ordinary (g) a-Gegenbauer] moments 
of the relevant DDs: 

= E j\ '^^ f2l, 2"-'= + x-^'^CiC^-'^ia) /^^(x, a; ^) da . (5.5) 



r(3/2)(n- 2fc)! 



Hence, each x"^Cf^^~^"^ (a) moment of f^^{x, a; fj.) is multiplicatively renormalizable and its evolution is governed 
by the anomalous dimension 7/+m (^i^- In Eq- ( ^.5| ), we took into account that a-DDs f{x,a) are always even in 
a, which gives an expansion of the Gegenbauer moments in powers of In the nonsinglet case, the Gegenbauer 
moments C„(^, /i) are nonzero for even n only. A similar represenation can be written for the Gegenbauer moments 
of the singlet quark distributions. In the latter case, the DD f^{x,a) is odd in x, and only odd Gegenbauer 
moments C,f (^, /i) do not vanish. 

Another simple case is the evolution of the gluon distributions in pure gluodynamics. Then the multiplicatively 
renormalizable operators with the same Lorentz spin n + 1 as in Eq. (tS.ll) are 



= z^^z-'izd+r-'G^^C'J^^^iz D /zd+) C . (5.6) 
Due to the symmetry properties of gluon DDs, only Gegenbauer moments 

C,?(C, /^) = £^ ct\{z/0 H'^iz, C; /i) dz (5.7) 
with odd n do not vanish. The Gegenbauer moment can also be written in terms of DDs: 

[(n-l)/2] I /-l-l^l rf Ohj_'l/0\ 

Two shifts: n ^ n — 1 and 3/2 5/2 in some sense compensate each other. Again, each combined x"^Cf^'^~^"^ (a) 
moment of f'~^(x,a) is multiplicatively renormalizable and its evolution is governed by the anomalous dimension 

Since the Gegenbauer polynomials Cf^^^™(a) are orthogonal with the weight (1 — q;^)™+^, evolution of the 
a;™-moments of DDs in both cases is given by the formula H 
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.1 oo 

U{a ;fi)= x'^fix, a;fi)dx^{l~ a^r+' ^ [log(A*/A)]-^'"+'/'''' , (5.9) 

•^-1 fc=0 

where the coefEcients Ami are proportional to x"^Cf^^~^"^{a) moments of DDs. A similar representation holds in 
the singlet case, with [log(/z/A)] "''''"+' ^'^^ substituted by a linear combination of terms involving [log(/i/ A)] "''''"+' 
and [log(/^/A)]-^"+'/^'' with singlet anomalous dimensions 7,^^; obtained by diagonalizing the coupled quark-gluon 
evolution equations 1^. 

The anomalous dimensions 7„ increase with raising n, and, hence, the mth x-moment of f{x, a; fj.) is asymp- 
totically dominated by the a-profile (1 — a^)™"*"^. Such a correlation between x- and a-dependences of /(x, a; 
is not something exotic. Take a DD which is constant in its support region. Then its a;'"-moment behaves like 
(1 — |a|)™~'"^, i.e., the width of the a profile decreases with increasing n. This result is easy to understand: due to 
the spectral condition \a\ < 1 — |a;|, the a;™ moments with larger m are dominated by regions which arc narrower 
in the a-direction. 

These observations suggests to try a model in which the moments /,„(q;;/i) has the asymptotic (1 — 0;^)™+^ 
profile even at non-asymptotic fj,. This is equivalent to assuming that all the combined moments x"^Cf^^^"^ (a) 
with I > vanish. Note that this assumption is stable wrt pQCD evolution. Since integrating [a ; fj,) over a one 
should get the moments fm{lJ^) of the forward density /(a;;^), the DD moments fm{.Oi]^) in this model are given 

by 

frnict;fl) = pm+lia)fmip) (5.10) 



where pm+i{a) is the normalized profile function (cf. Eq.(4^)). In explicit form: 



' x"^f{x,a;p)dx = -I^^±^Z^(1 - a2)'»+i f{z-p)z"^dz . (5.11) 
1 r(l/2) (m-l- 1)! 

In this relation, all the dependence on a can be trivially shifted to the Ihs of this equation, and we immediately 
see that /(x, a; p) in this model is a function of x/{l — a^): 

fix, a; n)^ F{x/(l - a^); p) e{0 < x/{l - a^) < 1) . (5.12) 

A direct relation between /(z, fj,) and F(u] fi) can be easily obtained using the basic fact that integrating f{x, a ; /i) 
over a one should get the forward density f{z,iJ,); e.g., for positive z we have 

f{z) = z f' —P^du. (5.13) 
This relation has the structure of the Abel equation. Solving it for F{u) we get 

^ Ju \JZ-U 

Thus, in this model, knowing the forward density /(z) one can calculate the double distribution function /(x, a) = 
F{xl(\-a^)). 

Note, however, that the model derived above violates the DD support condition -I- \a\ < 1: the restriction 
\x\ < 1 — defines a larger area. Hence, the model is only applicable in a situation when the difference between 
two spectral conditions can be neglected. A practically important case is the shape of H{x,£^) for small ^. Indeed, 
calculating H{x,£^) for small ^ one integrates the relevant DDs f{x) over practically vertical lines. If x is also 
small, both the correct \a\ < 1 — |a;| and model < 1 — |a;| conditions can be substituted by \a\ < 1. Now, if 
i 3> a slight deviation of the integration line from the vertical direction can be neglected and H(x,S.) can be 
approximated by the forward limit f{x). 

Specifying the ansatz for f{z), one can get an explicit expression for the model DD by calculating F{u) from Eq. 
( ^.14 ). However, in the simplest case when f{x) = Ax~'^ for small x, the result is evident without any calculation: 
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the DD /(x, a) which is a function of the ratio x/ (1 — a^) and reduces to Ax " after in integration over a must 
be given by 

J{x,a) = pa{a)f{x) 



where Pa{ct) is the normahzed profile function of Eq.(4.5) 



This DD is a particular case of the general factorized ansatz f{x,a) = pn{a)f{x) considered in the previous 
section. Its most nontrivial feature is the correlation n = a between the profile function parameter n and the power 
a characterizing the small- a; behavior of the forward distribution. 

Knowing DDs, the relevant SPDs H{x,£,) can be obtained in the standard way from f{x,a) for quarks and from 
xf'~^{x, a) in the case of gluons. In particular, the SPD enhancement factor Tl{() for small C, in this model is given 

by 



^c(0 r(2a + 2) 



fQ{C) r(a + 2)r(a + i) 

for quarks and by 

■Ff(C) _ r(2a + 2) 

C/G(C) "r(a + 3)r(a + i) 



(5.16) 



(5.17) 



for gluons. 

The use of the asymptotic profiles for DD moments /„(a) is the basic assumption of the model described above. 
However, if one is interested in SPDs for small ^, the impact of deviations of /„(«) from the asymptotic profile is 
suppressed. Even if the higher harmonics are present in /„(a), i.e., if the a;"~^'^C2^^^"^^'^(a) moments of f{x,a) 
are nonzero for k > 1 values, their contribution into the Gegenbauer moments C„(^,^) is strongly suppressed by 
^^'^ factors [see Eq.(^.5|)]. Hence, for small ^, the shape of H{x,£,) for a wide variety of model a-profiles is very 
close to that based on the asymptotic profile model. 

Absence of higher harmonics in /„(q;) is equivalent to absence of the ^-dependence in the Gegenbauer moments 
C„(^,/i). The assumption that C„(^,/i) do not depend on ^ was the starting point for the model of SPDs H{x,^) 
constructed in ref . |^ . Though the formalism of DDs was not used in ref. , both approaches lead to identical 
results: the final result of ||^ has the form of a DD representation for H{x,£,). In Appendix A, we also start with 
Cn{^,fJ-) = C„(0,^) and rederive the DD corresponding to the asymptotic profile model. 

VI. EVOLUTION ALGORITHM 

At one loop, evolution equations for nonforward parton distributions T^{x, p) can be written as 

dJ^^{x, fi) 



= Y,[wf{x,z):Fl{z,^ji^)dz (6.1) 



or in a "matrix" notation 



Using the explicit one-loop form of the effective coupling constant 
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477 

/3olog(MVA2) 



(6.3) 



{Po — 11 — 1^/) and the symbolic notations of (p^), one can present the formal solution for the set of evolution 
equations in the form of an expansion 



J^ciQ'') = cxp{L{Q^,Qo^) WA (».Fc(Qo') 



= E 



(6.4) 
(6.5) 



To get the n = 1 term SiTq = Wc_ ® TciQo'^) of this expansion, we evaluate numerically the convolution of the 
kernel with the initial distributions ^({Qo^)- To get the n = 2 term S2^c = (M^c)^ ® •^c('9o^), we convolute 
with the smoothly interpolated result of the first iteration SiJ-(^, and so on. After obtaining SiJ-(^, ^2^^^, etc., we 
construct evolved distributions !Fq{Q'^) for any desired value of Q^. Of course, the number of necessary iterations 
of with the initial distributions ^((Qo^) depends on the size of the expansion parameter 



i(Q',Qo') = -^iog 

Po 



log(QVA^) 
log(QoVA2) 



(6.6) 



When L is not very large, it is sufficient to calculate just one or two iterations. 



VII. EVOLUTION OF NONFORWARD QUARK DISTRIBUTIONS 

The evolution of skewed parton distributions was studied numerically in refs. |jl^,^,|2 14 1^. In this section, 
we perform the numerical evolution of SPDs using the algorithm described in the previous session and present the 



results illustrating evolution patterns for SPDs constructed using the factorizcd model (3.1) with different choices 
for the profile function h{x,a). 

As we discussed earlier, the use of the infinitely narrow profile function h(x, a) — 6{a) gives the simplest model in 
which OFPDs H{x,£^) coincide with forward distributions f{x). In terms of NFPDs, this modell looks less trivial. 
It gives F(°°)(a;,y) = S{y - x/2)f{x) for the y-DDs which results in (untilded) NFPDs given by shifted forward 
distributions 

For any monotonic function f{x) this gives NFPDs ^i;(X) which are larger in the region X > ( than their forward 
counterparts. Due to pQCD evolution, f{x,Q^) get steeper in small- a; region, i.e. the NFPDs become even more 
strongly enhanced. 

As noted in Section IV, the use of wide profile functions results in stronger enhancement for NFPDs in X > ( 
region. For Uustration see Figs.|[0. In Fig. ^we show C = 0.1 NFPDs obtained "by shift" from forward distributions 
f{x,Q^) taken at two values 1.5 GeV^ and 20 GeV'^. For comparison, we show also NFPDs obtained from 
f{x,Q^ = 1.5 GeV^) using the y-DDs with "asymptotic" profile functions Fj^'^(x, y) = 6y{l — x — y)f'^^{x)/{l — x)^ 
for nonsinglet quark distributions and F^^^^x, y) — 30y^(l — x — y)^/{l — x)^f'^{x) for gluon distibutions in pure 
gluodynamics. Specifically, we took /^^(a;,Q2 = 1.5 GeV^) = §x-i/2(i _ a.)3j^jjj f°{x,Q'^ = 1.5 GeV^) = 
0.4a;~°-^ln(l/x)(l - x)^ The NFPDS constructed in this way were then numerically evolved to — 20 GeV 
using the approach outlined in Section VI and kernels given Appendix B. The ratio of NFPDs obtained using these 
two models is shown in Fig.^. As expected, the ratios increase with Q^. 
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FIG. 6. YM gluon (left) and quark nonsinglet (right) nonforward distributions F^''^^{x) based on our model with has 
profiles versus "shifted" forward ones -(Yr^j2jf 1-^/2 ^ ' different scales; 1.5 GeV^ (dashed and solid curves) and 

20GeV^ (dotted and dash-dotted curves). 
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FIG. 7. Ratio of YM gluon (left) and nonsinglet quark (right) nonforward distributions F^'"'''{x) obtained from our 
model with has profiles to the "shifted" forward ones ■jj-r^j2)f ^ > different scales: 1.5 GeV^ (solid) and 

20GeV'^ (dashed); ( = 0.1. 

The same construction can be performed in the singlet case. Main observation here is large enhancement for 
singlet quark distributions at = 20GeV^, with the J^^(C)//"^(C/2) ratio being close to 1.8, see Fig.||. This is 
again in agreement with the estimates made in Section IV for awl. 
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FIG. 



Ratio of QCD gluon (left) and singlet quark (right) nonforward distributions ' (x) obtained from our model 
with has profiles to the "shifted" forward ones jj-r^jjjf ^fr^^i at two different scales: 1.5 GeV^ (solid) and 20 GeV^ 
(dashed); C = 0.1 - top, ( = 0.01 - middle, ( = 0.001 - bottom. 

In the above examples of nonsinglet quark distributions and gluon distributions in pure gluodynamics we took 
"asymptotic" profiles. It is interesting to test whether these profiles are really stable under pQCD evolution. To 
this end, we compared two models for = 20GeV^ distributions. First, we took the forward distribution evolved 
to — 20 GeV^ and constructed model NFPD using the "asymptotic" profile. Second way is to construct NFPD 
from = 1.5 GeV^ forward distribution using the asymptotic profile and then evolve NFPD to = 20GeV^ 
using nonforward kernels. Fig. ^ shows that the results obtained in the two ways are practically identical. 
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FIG. 9. Evolution of the YM gluon (left) and nonsinglet quark (right) nonforward distributions F^''^{x) obtained in our 
model with has profiles. Solid curves correspond to the initial distributions at = 1.5 GeV^. Dashed curves represent 
nonforward distributions evolved to — 20 GeV^, dash-dotted ones are obtained from the model with the same has profiles 
and with forward distributions evolved to = 20 GeV^ . 

However, if one takes profiles strongly differing from the "asymptotic" ones, the curves obtained in the two ways 
described above, visibly differ from each other, see Fig.|l0|. In the case of a wide profile function, the evolved NFPD 
looks like that constructed from evolved forward distribution but using a narrower profile. In other words, the 
pQCD evolution in this case narrows the profile function. Alternatively, if one starts with a too narrow profile, 
then the evolved NFPD resembles the model function constructed from evolved forward distribution but using a 
wider profile. 
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FIG. 10. The same as on previous set of figures (gluon distributions only) but with different profiles 

h{x,y) — N\{x){y{l ~ x ~ y))'^'- asymptotic (A = 2) - top, wide (A = 1) - middle, narrow (A = 6) - bottom. 

The study performed in Section V (see also Appendix A) demonstrated that at large there should be a 
correlation between the x-dependence of the forward distributions and the form of the profile function. Taking the 
GRV-type parametrization for gluon (with fg{x) ~ x~^-^ ln(l/a::)) and quark singlet distributions at — l.SGeV^, 
we again compared the — 20GeV^ curves constructed in two ways described above. A better agreement between 
the two models was observed for A = 1.5 rather than for A = 1.3. However, the A — 1.3 profile works perfectly if 
one takes the model with purely power-like behavior of the gluon distribution fg{x) ^ x~'^'^, see Fig.|lT|. 
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FIG. 11. Top and middle: the same as on previous set of figures but with other profiles h{x,y) = Nx{x){y{l — a; — y))'^- 
X = 1.3 - top, A = 1.5 - middle. Bottom figure represents the evolution of the model obtained with the profile parameter 
A = 1.3 and model forward distribution foi^) ~ l/a;" "^ (no logarithm). 



In the singlet case, it is more convenient to use "tilded" distributions defined on the [—1 + ( < X < 1] segment. 
These functions, shown in Fig.|l^, are symmetric or antisymmetric wrt the middle point X — C./2. 





-0 . 2 

FIG. 12. Evolution of the QCD gluon (left) and singlet quark (right) nonforward distributions F^'''{x) obtained our 
model with has profiles. Solid curves correspond to the initial distributions at = 1.5 GeV^. Dashed curves represent 
nonforward distributions evolved to = 20GeV^, dash-dotted ones are obtained from the model with the same has profiles 
and with forward distributions evolved to — 20 GeV^ . 



VIII. CONCLUSIONS 



In this paper, we discussed the structure of double and skewed distributions and their relation to the usual 
(forward) parton distributions. We emphasized that there are meson-exchange type terms in DDs and SPDs which 
are invisible in the forward limit. Theoretically, they can be modeled within approaches involving the mesons as 
elementary fields, e.g., this has been done in ref. Q within the chiral soliton model framework. There are also 
terms in SPDs and DDs which are related to usual distributions through reduction relations and which are, in this 
sense, visible from the forward limit. It is just these "forward visible" parts of DDs and SPDs which are the subject 
of our studies in the present paper. We proposed factorized models for double distributions in which one factor 
coincides with the usual (forward) parton distribution and another specifies the profile characterizing the spread of 
the longitudinal momentum transfer. By construction, these models satisfy reduction formulas. Then we used the 
factorized model DDs to construct skewed parton distributions and studied the skewedness dependence of SPDs. 
We established that, for small skewedness, the relevant parts of SPDs H{x,£,) can be obtained by averaging the 
usual parton densities /(i — ^a) with the weight p{a) coinciding with the profile function of the double distribution 
f{x, a) at small x. This result allows one to get estimates for the ratio of SPDs taken at the border point x = ^ 
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(or X = C) and usual parton densities taken at x = xsj- This ratio is larger than 1 in all cases which we have 
considered, i.e., SPDs in this sense are enhanced compared to the forward densities. We found that, for small xbj, 
the enhancement factor is directly related to the parameter a characterizing the effective power behavior a;~" of 
the usual parton densities. We gave an explicit expression for the enhancement factor which involves a and the 
parameter b characterizing the effective power profile (1 — a^)^ of the relevant DD. Furthermore, we demonstrated 
that if the moments /„(a) of the a-DDs have the asymptotic (1 — q;^)"+^ profile, then the parameters a and 
b are correlated: b = a, i.e., the a-profile of f{x,a) for small x is completely determined in this case by the 
small- a; behavior of the usual parton distribution. We also demonstrated that in the case of small skewedness ^, 
the deviation of fn{o:) from its asymptotic form leads to very small O(C^) corrections only. This means that the 
"forward visible" parts of SPDs can be rather accurately deduced from the purely forward usual parton densities. 
The absence of higher harmonics in fnio:) is equivalent to the absence of ^-dependence in the Gegenbauer moments 
of SPDs, which is the starting point of the model for SPDs constructed in the recent paper by Shuvaev et al. For 
this reason, the results based on the asymptotic profile model (APM) developed in the present paper coincide with 
those of ref. though its authors refrained from using DDs in their approach. Finally, we performed a numerical 
investigation of the evolution patterns of SPDs and gave interpretation of the results of these studies within the 
DD framework. In particular, if the initial profile of DDs is too wide, the pQCD evolution makes it more narrow 
and the profile widens if initially it was too narrow. Our numerical results also support the expectation that if one 
takes SPD derived from the forward distribution by APM prescription and evolves it using nonforward evolution 
kernels, the result coincides with SPD which is APM-derived from evolved forward distribution. Summarizing, 
our results based on numerical evolution of SPDs are in full accord with theoretical expectations based on DD 
formalism. The next step in this investigation should be a direct numerical evolution of double distributions. 
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APPENDIX A: ASYMPTOTIC PROFILE MODEL AND ^-INDEPENDENT GEGENBAUER 

MOMENTS 

The lowest k = term of the ^■^'^ expansion for the Gegenbauer moments of the SPDs involves the DDs f{x,a) 
integrated over a with Cq^^^^ (a) — 1, i.e. the ordinary (forward) parton distributions: 

C«({ = 0,rt = 2»S|i5Za |'^/0(,;p),-<ix. (Al) 

- "> - 2" mii^. L ■ 

The higher terms of the ^^'^ expansion are small for small ^. The approximation in which the k terms of 
the ^^'^ expansion for the Gegenbauer moments are neglected (i.e., Gegenbauer moments C„(^,/i) are treated as ^- 
independent) can be translated into a model for the double distribution /(x, a). The advantage of this model is that 
its basic assumption is stable with respect to evolution: if the Gegenbauer moments C„(^,/io) a-re ^-independent 
at some low normalization point ^q, this property is preserved for any higher /i. The only impact of the pQCD 
evolution is the change of the relevant forward parton distribution from f{x,fio) to f{x,fi). 

To derive the explicit expression for f(x, a) in this model, we will use the expansion of the light-cone operator 
'ijj{—z/2)X°'zip{z/2) over the multiplicatively renormalizable operators 0„ (see [|l8| ) 
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V;(-z/2)A»£V(^/2) = E(-l)" / (1 - a^T+^0,Aaz/2)da . (A3) 

n=0 * 

Inserting it into the nonforward matrix element, we obtain 
(P - r/2 I V;(-^/2)A'^z^(z/2) | P + r/2 ) = j' ^ f;(-l)"-ij^±^(l _ a2)"+i(-zP.)"C„(^, /.) da • 



(A4) 



If the Gegenbauer moments C„(^,/i) are approximated by their ^ = values, Eq.(A4) can be transformed into 
the representation of the nonforward matrix element in terms of double distributions. Namely, the nth moment 



»i 

/„(a;/x)= / x"-f{x,a;fi)dx (A5) 



of the a-DD f{x,a; fi) is then given by 

fnia ; A*) = r(V2H»^+V ^ ~ /'i ^^^"'^^ ■ ^^^^ 

The factor relating /n(a; /i) and fnifJ-) is just the normalized profile function pn+i{a) (see Section ?). In the 
nonsinglet case only moments with even n are nonzero, while in the singlet case only those with odd n do not 
vanish. Below, we will construct the x > parts of DDs, the trivial (anti)symmetrization can be performed at the 
end. 

Incorporating the inverse Mellin transformation, one can obtain the kernel K{x, a; z) which relates DDs /(x, a) 
in this model with the usual forward distributions f{z). However, inverting the representation ( |A6| ) "as is", one 
would get the expression 

1 

f{x,a)^-^ f [z{l-a')/x-l]-'/'f{z)-, (A7) 
2tt J z 

whose rhs has a suspicious overall sign. Furthermore, the integral over z diverges at the end-point z — x/ (1 — a^). 
These inconsistencies indicate that the implied interchange of the inverse Mellin transformation and the z-integral 
is not justified. To get a less singular kernel, one can try to add an 0(l/n) factor in the expression for the moments 
/„(a; /u), e.g., convert l/(n + 1)! into l/(n + 2)! (after such a change, the inverse Mellin transform is still easily 
doable analytically). To this end, we use the relation 







/(z; ^i)z"dz = -—^ I z"+^ — ( ] dz (A8) 



which holds if (a) the function f{z; ji) vanishes aX z — \ (this is always true) and {b) z"^^^ f{z) vanishes at z = 0. 
The latter requirement is evidently satisfied in the nonsinglet case for all n > 0. The singlet distributions are more 
singular, but we need only their n > 1 moments, i.e., the restriction (6) is satisfied again. Using Eq. ( |A8| ) produces 
the kernel which connects f{x,a) with [f{z)/z]': 



1 

/(^,«)-- ^(^!^2) / Hl-a')/x-ir'^'[f{z)/z]' dz . (A9) 

x/{l-a^) 



This result coincides with Eq.(5.14). The spectral condition x/z < 1 — relating the "original" fraction z and 
the "produced" fraction x is analogous to the momentum ordering x < z in the DGLAP equation: the produced 
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fraction cannot be larger than the original one. In the present case, if the parton also takes some nonzero fraction 
a of the momentum transfer, the allowed values of x cannot exceed z{l — a^). 
For gluons, we combine the expansion for the bilocal operator 



3(2n + 5) 



n=0 



22»+3(n + 2)! 



{1 ~ a^)''+^0^+i{az/2) da 



(AlO) 



and the expression (|A^ ) for the Gegenbauer moments at ^ = 0. The resulting relation for the moments /„ (Q;;/i) 



has the form identical to Eq.(A6) derived in the quark case. 

APPENDIX B: EVOLUTION EQUATIONS 



In Ref. 1^, the kernels W^''{x, z) were obtained via 



B^^u, v)6 {x~uz + v{z - 0) e{u + v < 1) du dv 



from the evolution kernels for the light-ray operators fl^,^ : 

B11 = ^Cf ( 1 + Siu)[v/v]+ + S{v)\u/u]+ - l5{u)6{v) 
n \ 2 

B31 = (2 + 5{u)5{v)) 

TT 
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'-Nf{l + Auv ~u~v) 



f3o 



—Nc 4(1 + 3uv-u-v) + 777r'^(")^M + ^i^)^ 







'1' 












V 


+ 


u 





(Bl) 

(B2) 

(B3) 
(B4) 

(B5) 



In the nonsinglet case, only the qq kernel is needed. Using Eqs. (B2),(B1), one can easily derive the rules allowing 
to transform each of the four terms contained inside the bracket in Eq. (B2) 



e{x<0 J ^-^^Tl{z)j+ J !^9ix<C)~ + eix>C)^\T^,{z)^d^, 

X 
miii{a: / Ct^/C} 



Siu) 



dv 



(x<C) 



+ 0{x > C) 




q I X-VQ 



TKx) 



5{v) 



'W 


.-I 


.u. 










du 



dv 



dv 



[l + log(min{x/C,x/C})]^^'(x) 



+ [l + log(x/C)]^^^(x) 



-+[l + \og{x/0]mx)} ^d2 



Tl{x/u) ~ Tlix) — + (1 + log(a;)) Tl{x) = d^, 



-^S{u)Siv)^-^:F^^{x)^d,. 



(B6) 



(B7) 

(B8) 
(B9) 



In terms of these contributions, the total qq part of (6.1) is given by 
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® PUll) = —Cpidl +d2 + d3 + d4). 



(BIO) 



Expressions for di and c?4 were obtained directly by performing integrations over u, v of the terms in Eq. (B2) 
with 1 and {—^S{u)S{v)). To obtain expressions for ^2 and ds, we changed the order of integrations. For example, 
to obtain da, we first take the integral over v and then over z: 



\Jo Jo 



\"'0 



S{x — uz + v{z — ^))(5(u) — ^(;(^) du dv j dz 



1 du 



S{x — uz) — !F(^{z) du \ dz — / ^9{x/u<l) — !Fq{x/u) 



du 



J-q{x/u) ^c(^) 



■ du 



du r 



T'^{x/u) - T'^ix) + (1 + log(5)) T'^ix). 



(Bll) 



The resulting form for integrals in d2 and ds is particularly convenient for numerical calculations, because the 
integrand in Eqs. (B7), (B8) is explicitly regular as u,v ^ 0. 
In a similar way, we represent the gg -kernel part as 



W^f ^ THp^) = —N.idl + dl + d',+ dl). 



(B12) 




(C(2 + x- z)-x^ - z{2~ x)) T^iz) dz + 



2x C — X 



dz 

z 



Tl{z) 



dz 

z 



(B13) 



-TUx) 



Tlix) 



dv 

V 



-+\og{xlOmx) 



dv 



\og{x/OTl(x)) , 



4- 



{uTl{x/u) - Tl{x)) ^ + \Qg{x)Tl{x), 



(B14) 



(B15) 



(B16) 



The qg and gq parts of the evolution equations (3.2) for untilded NFPDs in the region where z < C cannot be 
unambiguously reconstructed from the light-ray kernels. We use the form suggested in [H, which leads to 



{x > C) 



C — X f Z ~ X 



C J z-C 



x{x - C) 
z{z - C) 



dz 

z 



+ 0{x < C) 



iC^xy 



1 



2z x(C^ 



JF?(z) dz 



(B17) 
(B18) 
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1 

-, 2z 2(C-a;) 



^2 ; ^ l + j+^^^^^^]^!{z)dz}, (B19) 



(C-x) /• 1 / 2x x(C-z) 



(B20) 



+ eix<o{^Y^ y ___|^i__„4^^)^^^(z)dz (B21) 



1 

2x z{C-x) 



X r 1 



C C 



In the matrix notation 

W'" ® - —Cfs^i, (B23) 

W-«s ® ^^3(^2) = ^TVfs'Js. (B24) 
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